Indirect Evidence for Levy Walks in Squeeze Film Damping 



O. 






cr 
i , 



> 

vo 

On 

in 

o 
o 



X 



S. Schlamminger, C. A. Hagedorn, J. H. Gundlach 

Center for Experimental Nuclear and Particle Astrophysics, 

University of Washington, Seattle, WA 98195 USA 

Molecular flow gas damping of mechanical motion in confined geometries, and its associated noise, 
is important in a variety of fields, including precision measurement, gravitational wave detection, 
and MEMS devices. We used two torsion balance instruments to measure the strength and distance- 
dependence of 'squeeze film' damping. Measured quality factors derived from free decay of oscillation 
are consistent with gas particle superdiffusion in Levy walks and inconsistent with those expected 
from traditional Gaussian random walk particle motion. The distance-dependence of squeeze film 
damping observed in our experiments is in agreement with a parameter-free Monte Carlo simulation. 
The squeeze film damping of the motion of a plate suspended a distance d away from a parallel surface 
scales with a fractional power between d~ l and d~ . 

PACS numbers: 05.40.-a,07.10.Pz, 95.55.Ym,07.30.-t 



The damping of mechanical motion by residual gas par- 
ticles in vacuum has been the subject of measurement, 
discussion, and controversy [l|-l5j]. 'Squeeze film damping' 
is enhanced gas damping that occurs due to molecular 
flow in constrained geometries [6-9]. Quantitative un- 
derstanding of this damping is essential to a wide range 
of subjects: The quality factor of many MEMS oscilla- 
tors is limited by squeeze film effects. In future grav- 
itational wave interferometers, test masses will be sus- 
pended within millimeters of another surface, constrain- 
ing gas flow. The fluctuation-dissipation theorem dic- 
tates that squeeze film damping generates noise that may 
limit these experiments [10] . 

Noteworthy are two recent publications on squeeze film 
damping: In Ref. [fj the damping of three topologically 
different MEMS oscillators suspended a fixed distance 
from a substrate was measured as function of residual gas 
pressure. Ref. [7] also measured the damping as function 
of gas pressure, using two torsion balances centered in 
closely-spaced housings. Both groups measured greater 
damping than expected from free gas damping, and de- 
veloped Monte Carlo simulations that agreed with their 
experimental data. Ref. [7] drew attention to the impor- 
tance of constrained gas damping to gravitational wave 
detectors. As their geometry models the planned LISA 
gravitational wave detector, which is sensitive to both 
transverse and longitudinal gas flow, they did not repli- 
cate the planar geometry of existing and planned ground- 
based interferometers. These experiments and the failure 
of analytic theory to match their results motivated our 
measurement of the distance dependence of squeeze film 
damping with a simple planar geometry. 

We used two different torsion balances to investigate 
mechanical damping in rarefied gas. We provide experi- 
mental data characterizing the distance and pressure de- 
pendence of squeeze film damping using the simple and 
well controlled geometries of two adjacent surfaces. We 
compare our data to analytic theory and to Monte Carlo 
simulations. The simulations indicate that Gaussian ran- 
dom walks incorrectly describe the diffusion of gas par- 
ticles in narrow gaps. Instead, squeeze film damping, 
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Figure 1. Schematic drawing of the LTA torsion pendulum. 
The copper plate is mounted on a translation stage allowing 
gaps up to 9 mm. The pendulum is suspended by a quartz 
fiber with Q > 50000. 



including distance-dependence, is better described by su- 
perdiffusive particle transport in Levy walks. 

Each of our two torsion balance instruments consists 
primarily of a thin vertical plate pendulum suspended 
from a torsion fiber. Each pendulum is parallel to a sta- 
tionary surface situated a distance d away. We deter- 
mined the damping by measuring the mechanical quality 
factor, Q, of each pendulum's free oscillation as a func- 
tion of d and as a function of the residual gas pressure 
P. 

The LISA Test Apparatus (LTA) torsion balance [HI 
depicted in Fig. [1] was constructed specifically to investi- 
gate noise sources in LISA. The plate of the LTA torsion 
pendulum is a 0.45 mm thick Si wafer, 114.4 mm wide 
and 38.1 mm tall. The plate has a small riser section con- 
necting it to a perpendicular Al bar designed to minimize 
its sensitivity to fluctuating gravitational fields caused by 
nearby human activity. The pendulum assembly is sus- 
pended from a quartz fiber with intrinsic Q > 50000. 
The assembly has a free period of 120.89 s and a mo- 
ment of inertia of 9.7 x 10 -6 kgm 2 . The pendulum plate 
is positioned parallel to a thick Cu plate. The Cu plate is 
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Figure 2. Schematic drawing of the NTA torsion pendulum. 
The pendulum is suspended by a tungsten fiber with intrinsic 
Q ~ 4000. The pendulum can be translated to change the 
gap between pendulum and foil. 



divided into two halves; both are electrically grounded. 
The pendulum to Cu-plate gap, d, can be varied between 
and 9 mm. The Cu plate and pendulum are sputter- 
coated with Au to reduce contact potentials and patch 
effects [12]. The total electrical charge on the pendulum 
was reduced below 0.5 pC [l3j; our results were not af- 
fected by deliberate increases in charge. The torsion bal- 
ance is housed within a stainless-steel vacuum chamber. 
The vacuum is maintained by a small turbo-molecular 
pump with an extra molecular drag stage backed by a 
small diaphragm pump. 

Our second instrument, the 'Newton Test Apparatus' 
(NTA), is depicted in Fig. H A 2.00 mm thick, 42.18 
mm wide, and 30.81 mm tall metal plate is suspended 
from a tungsten fiber near a tightly-stretched 12 /jm thick 
flat beryllium-copper foil. The foil-pendulum gap, d, is 
varied by moving the torsion fiber suspension point. The 
pendulum has a free period of 85 s and a moment of 
inertia of 3.1 x 10~ 6 kgm 2 . Like in the LTA, both foil 
and pendulum are coated with Au. The torsion balance 
is housed within a stainless-steel vacuum bell jar and the 
vacuum is maintained by a turbo-molecular pump. 

The torsional motion of each pendulum is monitored 
with a laser autocollimator. 

To determine the size of the gap and to assure proper 
alignment between the pendulum plates and their adja- 
cent surfaces in LTA and NTA, several methods were em- 
ployed: (1) Calibrated photographs were taken through 
vacuum windows. (2) Each pendulum was allowed to 
twist in both directions until the pendulum's outer edge 
touched the adjacent surface. The gap was determined 
by the maximum twist angles and the size of the pendu- 
lum plate. (3) In the LTA, the capacitance between the 
two halves of the copper plate was measured as a function 
of d and was fit to a capacitance model. Discrepancies 
between the distance measurements determined the un- 
certainty ad = 0.1 mm in LTA, 15 fira in NTA. The 
parallelism, about the horizontal axis, of the pendulums 



w.r.t the adjacent surfaces was better than 0.5 mrad in 
LTA and 30 /irad in NTA. The parallelism about the ver- 
tical axis was better than 0.15 mrad in LTA and 0.6 mrad 
in NTA. Systematic and flatness distance uncertainties 
were < 45 /xm in LTA and < 23 fj,m in NTA. 

To improve confidence in our pressure measurement, 
we connected five ion gauges from three different man- 
ufacturers and a variable leak valve to a single small 
vacuum chamber. With the exception of one suspect 
gauge, the gauges agreed to within 30% over the range 
of 10~ 4 -2.6 x 10~ 2 Pa. From this set, a cold-cathode 
gauge was used for all of the LTA measurements and a 
hot cathode gauge was used for all of the NTA measure- 
ments. We calibrated the hot-cathode ion gauge with a 
simple pressure-division apparatus. With it, we verified 
the gauge linearity and calibrated the gauge absolutely 
to within 2 x 10~ 3 Pa. Our experiments were carried out 
at a temperature of 297 ± 1 K. 

In the LTA, the pressure was adjusted by changing the 
speed of the turbo pump and allowing the apparatus to 
come into equilibrium. Data were taken at three different 
pressures. The pressure at the gauge was ~ 10% smaller 
than the pressure in the chamber due to gas pumping 
impedances. We measured the composition of the gas 
in the chamber with a residual gas analyzer at all three 
working pressures. H2O (18 amu) and N2 (28 amu) to- 
gether accounted for at least 85% of the residual gas. 
At the working pressures of 2 x 10~ 3 , 1 x 10~ 4 , and 
4 x 10~ 5 Pa, the partial pressure ratios (pis'P28) were 
0.4:1, 2:1, and 8:1, respectively. Nitrogen is pumped more 
efficiently at high turbo pump speeds. 

In the NTA, measurements were made by pumping 
the vacuum vessel below 10~ 4 Pa, closing the turbo 
pump's roughing valve, and turning the turbo pump 
off. Q measurements began after the turbo had stopped 
(P > 9 x 10 -4 Pa). The pressure in the chamber rose 
at a constant rate, most likely due to outgassing from 
the chamber walls. Data were taken with pressure rises 
of 1.3-7 x 10~ 7 Pa/s until the pressure reached 1.3- 
2 x 10~ 2 Pa. We assume that the residual gas was pri- 
marily particles with mass 18 amu, as we expect the out- 
gassing to be dominated by water. 

We measured the free decay of each pendulum's tor- 
sional mode to determine Q. The LTA torsion pendu- 
lum was excited to an oscillation amplitude of 0.9 mrad. 
The pendulum motion was recorded for > 10000 s. 
The data were cut into segments two oscillation peri- 
ods in length. The data in each cut were fit to 9{t) = 
A cos(2ir ft) + B sm(2ir ft) + C. The extracted amplitudes 
were then fit to an exponential decay to extract Q, where 
the amplitude 6(t) = o e _7r/ * /Q = vW+W ■ 
Errors of 20% were deemed to be conservative overes- 
timates in all cases. In the NTA, data recordings of 
>1000 s, with initial amplitude 0.2-4 mrad, were fit to 
9(t) = e~* t/TQ (A cos(27rf/T) + B sin(27rf/T)) + C + Dt, 
where A, B, C, D, T, and Q were fit simultaneously. This 
procedure allowed measurements while the gas pressure 
was slowly varied. Seismic noise and other excitations 



conservatively limit our ability to measure Q to < 50000 
in the LTA and < 1500 in the NTA. 

The inverse of the measured quality factors, as func- 
tions of pressure and pendulum-plate separation, are pre- 
sented in Figs. [4] and [5] 

In the molecular flow regime, where the particle mean 
free path is much larger than the gap and particle interac- 
tions may be neglected, gas damping is primarily due to 
'unconstrained gas damping' and 'squeeze film damping'. 
All of our data were taken in this regime. 

In unconstrained gas damping, the motion of an ob- 
ject is retarded by collisions with ambient gas particles. 
A rectangular plate of dimensions a x b moving perpen- 
dicular to its surface with velocity u through a dilute 
gas is retarded by the force: F(u) = —aPab(v)~ l u, with 
(v) = y'WtBT/Jjrm) , m the gas particle mass, the ambi- 
ent pressure P and temperature T. Three values for a 
can be found in the literature: a — 16/7T 0,@|, a = 16 Q, 
and a = 24 [l|. Unconstrained gas damping depends lin- 
early on the gas pressure but is independent of the gap 
to nearby surfaces. 

Squeeze film damping of two adjacent plates moving 
towards each other arises because the gas between the 
plates is temporarily compressed until the "trapped" gas 
can diffuse out of the gap and restore equilibrium with 
the surrounding gas. To compute the dissipative force on 
the plate we start with Suijlen's derivation J6J. Periodic 
motion of the plate with amplitude zq <C d and angular- 
frequency u <C (v)/d will cause the distance-dependent 
squeeze film force on the moving plate: 
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where the diffusion time, r, for particles to enter or leave 
the gap is discussed in detail below. 

For a torsional oscillator, the damping forces must be 
converted to torques. We approximate the torsion pen- 
dulum as two a x b plates with lever arm a/2. We expect 
that a more precise model would give only small correc- 
tions to this approximation. For a torsion balance with 
resonant frequency /q <C 1/t and nronrent of inertia /, 
we find that the inverse of the quality factor Q is: 
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Since Suijlen et al. [6] and Cavalleri et al. \7\ found 
that their experimental results and simulation disagreed 
with Gaussian theoretical expectations, we undertook 
our own simulation effort to deternrine the diffusion tinre 
r. Wc simulated gas particles traveling interaction- 
free in the gap. Initially, particles were randomly 
distributed throughout the gap volume and assigned 
three-dimensional velocities according to the Maxwell- 
Boltzmann distribution. Upon collision with a wall, the 
particles were emitted with a polar angular distribution 
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Figure 3. Simulated trajectory of a gas particle in the NTA 
experiment for a 0.1 mm gap. The occasional large steps are 
a consequence of the infinite variance of the step-size distri- 
bution. These steps, which lead to superdiffusion, and scale- 
invariance are characteristic of Levy walks. 



p{9) and a Maxwellian velocity distribution. The angular 
and velocity distributions of desorbed particles are poorly 
understood. A review, [14,], suggests that p{9) oc cos™ 9 
with n > 1, where the traditional choice, n = 1, is mo- 
tivated by time-reversal symmetry. The particles were 
tracked until they left the volume, and the time at which 
they left the gap was recorded. The average of these 
times, for 10000 trials, was taken as the simulated diffu- 
sion time. For the setup for which Suijlen et al. provide a 
complete description of their geonretry, our results agree. 

Simulated particle tracks (Fig. [3]) show occasional very 
long steps. These long steps are essential for understand- 
ing squeeze film damping. 

Assuming thermal equilibrium and p(9) oc cos 9, the 
probability for a particle to make a step of length r s with 
a duration t s is given by 



Pi(r Sl t s )dr s dt s 
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The length of the step, r s , is measured parallel to the 
plates. r s and t s are not independent, but the distri- 
bution of each can be obtained by integrating pi(r s ,t s ) 
over the other. For example, the steps are distributed 
as pi(r s )dr = 2d 2 r s /(d 2 + r 2 ) 2 dr s , consistent with [6]. 
The step-size and time-of-flight distributions have infi- 
nite variance, indicative of Levy flights/walks. In a Levy 
flight the tinre to make a step is independent of step 
length. In a Levy walk the particle travels with finite 
velocity [15H18J . as in a physical gas. The defining fea- 
ture of Levy flights/walks is that repeated application of 
the step-size distribution leads to a stable, scale-invariant 
distribution with infinite variance. The Lorentzian is an 
example of such a distribution. Our Monte Carlo sim- 
ulations show that pn{t) becomes scale invariant after 
a few bounces. Levy walks can lead to superdiffusive 
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Figure 4. The dissipation factor (1/Q) as a function of plate- 
pendulum separation at constant pressure in the LTA. The 
lines are calculated using Eq.[5]with the diffusion time r deter- 
mined by the parameter-free Monte Carlo simulation (solid) 
and the Gaussian diffusion equation (dashed). The curves 
flatten with increasing gap as unconstrained gas damping be- 
gins to dominate. The residual gas composition changes with 
pressure; the curves are calculated using the five measured 
leading gas species at each pressure. 



behavior [Fa , H8j |, i.e. the mean square displacement of 
an ensemble of particles released at r = increases as 
(r 2 N ) oc V , with 7 > 1. Our Monte Carlo simulation 
of gas particle trajectories trapped between two infinite 
planes shows 7 « 1.15, independent of d, for t 3> d/(v), 
i.e. after a few bounces. (r%), in contrast to (r^), re- 
mains finite because (rj^) is evaluated using the position 
of all particles in free flight between collisions with the 
two planes. 

The smooth curves in figures U and [3] are computed us- 
ing Eq. [21 with a = 16. The dashed lines are calculated 
with the Gaussian diffusion time, r = 8ab/(ir 3 d(v)) Q, 
and the solid lines are computed using our Monte Carlo 
diffusion times. The dotted line in Fig. [5] is computed 
using the same simulation but with an angular distribu- 
tion of emitted particles p(9) oc cos 2 9. As the damping 
is expected to be linear in pressure, we also plot 1/(QP) 
for the NTA. We cannot do so for the LTA results, as 
the residual gas species change with pressure. The lines 
in Fig. [4] are calculated with the measured residual gas 
composition. 

Our squeeze film measurements are inconsistent with 
Gaussian diffusion. In Eq. [3J for small d, Q is propor- 
tional to cr, where j3 depends on the diffusion time, r(d). 
For r oc d~ x (Gaussian diffusion), j3 — 2 and for constant 
r, $ = 1. We find (3 = 1.4 ± 0.3 in LTA and 1.6 ± 0.3 
in NTA. Our Monte Carlo simulations predict j3 = 1.73 
in LTA and 1.85 in NTA. The data are consistent with 
our simulations using p{9) oc cos™ 9 for n = 1 and clearly 
inconsistent for n > 2. The measured Q is inversely pro- 
portional to gas pressure, as predicted in Eq. [5] 

We have made the first rigorous measurement of the 
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Figure 5. Dissipation factor as a function of pendulum-foil 
separation in the NTA. The lines are calculated as in Fig. [?] 
The lower panel shows the dissipation factor divided by the 
pressure. These data agree well with our Monte Carlo predic- 
tion (solid line), which has superdiffusive Levy-like behavior, 
The data are inconsistent with Gaussian theory (dashed line) 
and with a desorption angular distribution p(6) oc cos n 9, for 
n > 2 (dotted line). 



distance dependence of squeeze film damping by adjust- 
ing the gap between an oscillator and a flat plane. Exper- 
iments in two separate instruments were performed using 
planar geometries such that the gas flow is restricted by a 
single gap. This geometry is analogous to the geometry 
in existing gravitational wave interferometers and sim- 
pler than the geometry of the experiments decribed in 
Ref. [7| • We point to the connection between squeeze film 
damping and the theory of Levy walks. The infinite vari- 
ance of the step-size distribution of Levy walks leads to 
superdiffusion, which is consistent with our data. These 
measurements are simple examples of microscopic Levy 
phenomena in nature. In addition, our results support 
a cosine angular distribution of re-emission of adsorbed 
gas particles. As predicted by our Monte Carlo simula- 
tion, the measured squeeze film damping is smaller and 
less dependent on gap size than expected by Gaussian 
theory. The comprehensive data set presented here can 
be used to verify any Monte Carlo simulation of squeeze 
film damping. Squeeze film damping must be considered 
as a limitation of next-generation gravitational wave de- 
tectors. The design of these detectors requires precise 
knowledge of the strength and distance dependence of 
squeeze film damping. 
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